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ON THE MULTILINEAR HAUSDORFF PROBLEM OF 

MOMENTS 

A. IBORT AND P. LINARES, J.G. LLAVONA 

Abstract. Given a multi-index sequence /ik, k = (k\, . . . , fc„) € 
Nq, necessary and sufficient conditions are given for the existence 
f- ^ , of a regular Borel polymeasure 7 on the unit interval / = [0, 1] 

£NJ ' such that ^ k = J* „ t^ 1 ® ■ • • ® £*" 7. This problem will be called 

^ . the weak multilinear Hausdorff problem of moments for fi^. Com- 

parison with classical results will allow us to relate the weak mul- 
tilinear Hausdorff problem with the multivariate Hausdorff pro- 
blem. A solution to the strong multilinear Hausdorff problem of 
moments will be provided by exhibiting necessary and sufficient 
conditions for the existence of a Radon measure /1 on [0, 1] such 
that Lp(fi,...,f n ) = f T fi(t) ••■f n (t) V>{dt) where i M is the n- 
~^<i linear moment functional on the space of continuous functions on 

I^H . the unit interval defined by the sequence /ik- Finally the previ- 

ew ' ous results will be used to provide a characterization of a class 

of weakly harmonizable stochastic processes with bimeasures sup- 
ported on compact sets. 
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\o 

The moment functional L M associated to a sequence /j,k, k G No of 
real numbers is the element in the (algebraic) dual of the space of poly- 
nomials R[t] defined by L^(p) = J2 k > Pk^k, where pit) = ^k>oPk tk £ 
M.[t] is an arbitrary polynomial. Given an interval /cl the classical 
problem of moments for the sequence fi k asks for the integrality of the 
^ linear operator L M , that is, under what conditions there exists a (pos- 

sibly signed) Radon measure //on/cK such that L^(t k ) = fjt k dfj,{t), 
k = 0, 1 ... . If / is the unit interval [0, 1] the problem of moments is 
known as the Hausdorff moment problem. 

The well-known solution to the classical Hausdorff problem (see for 
instance |Sh70] and references therein) establishes that such a measure 
\x exists provided that there is a constant C such that: 

k 

(I- 1 ) 2-^ \\k;m)\ < C, 

m=0 



a 



1. Introduction 



Key words and phrases. Problem of moments, multilinear, polymeasures, second 
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for all k = 0, 1, . . ., where \ {k , m) = ( ^ J L M (£ m (l - i) fc ~ m ). 

There are various natural extensions of the moment problem to the 
multilinear case. Given the multi-index sequence fi^, k = (hi, ... , k n ) e 
Nq, we will denote as before by L^ the n-linear functional defined on 
the space of polynomials M[t] by: 

(1.2) L,(t k \...,t k -)=fi k Vk=(k 1 ,...,k n )eN l . 

Thus the strong multilinear Hausdorff problem of moments will consist 
in determining under what conditions there exists a (possibly signed) 
Radon measure \i on [0, 1] such that 
(1.3) 

L^ Pl (t),..., Pn (t)) = / Pl {t)--- Pn {t)n{dt), \/ Pl (t),..., Pn (t)eR[t]. 

Jo 

A weaker version of the multilinear Hausdorff moment problem, the 

classical (multivariate) Hausdorff moment problem, can be stated by 

demanding the existence of a (possibly signed) Radon measure fi on 

[0, l] n C R n such that: 



(1.4) L ll (p 1 (t 1 ),...,p n (t n ))= Pi(t 1 )"-p n (t n )fi(dti,...,dt n ). 

J i" 

Integrality properties of bilinear functionals related to the notion 
of the total variation of functions in many variables were studied by 
Morse and Transue |Mo49j . These questions led to the concept of C- 
bimeasures |Mo56j . These results were also deeply rooted in the pro- 
blem of studying the structure of stationary stochastic processess (see 
for instance |Ra82] and references therein). The notion of bimeasures, 
and polymeasures in general, provide a natural framework to answer 
these questions. Thus we are naturally led to consider an even weaker 
version of the classical multilinear Hausdorff problem of moments. We 
will say that /ik satisfies the weak multilinear Hausdorff problem of 
moments if there exists a polymeasure 7 on Bo[0, 1] x ■ ■ • x Bo[0, 1] 
such that: 

(1.5) /i k = / t?®---®t£>y(dt 1 ,...,dt n ), Vk=(k u ...,k n )eNl 

As it turns out, the solution to the weak multilinear Hausdorff pro- 
blem of moments is given by a nontrivial generalization of condition 
eq. ( 11.1)1 as it will be proved in section 2, Thm. 12.21 This condition is 
different from the characterization obtained in the analogous weak mul- 
tilinear trigonometric problem of moments [IblOj . Using these ideas the 
classical condition ( II. ip and the solution to the classical multivariate 
Hausdorff problem are easy consequences of the general properties of 
polymeasures on compact sets as it will be discussed briefly at the end 
of section 2. The strong multilinear Hausdorff problem of moments will 
be solved in section 3 by using recent results on integral representations 
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of orthogonally additive polynomials on Banach lattices. In addition 
the existence of the integrating measure will be completely character- 
ized as an algebraic property of multimoment sequences, Thm. 13.11 
Finally, as an application of the previous results we will provide a new 
characterization of a class of weakly harmonizable stochastic processes 
in section HI 

2. The weak multilinear Haussdorff moment problem and 
a multilinear Riesz theorem 

A polymeasure 7 on the cr-algebras Ei,...,S n is a separately a- 
additive function on the cartesian product of Ei,...,E n [Do87j (we 
will consider here only real or complex polymeasures). The variation 
of the polymeasure 7 is defined as the set function v(j): Si x • • • x E n — > 

[0,+oo]: 

v( 1 )(A l ,...,A n )^su V lj2---J^\ 1 (Al\...,A k n -) 

where the supremum is taken over all finite partitions {A l l } r k l =1 of the 
set A\ G E/. The semivariation ||7|| : Ei x • • • x E n — > [0, +00] of the 
polymeasure 7 is defined as: 
(2.1) 

{ Ji r n 

where the supremum is taken over all finite partitions {A l l } r k l =1 of the 

set Ai G E/ and all collections of numbers {a l l } r k l =1 such that \a h l \ < 1. 
In the linear case n = 1 the semivariation and variation of a measure 
coincide. 

An integral denoted as f f\ <g) . . . ® f n 7, can be constructed for poly- 
measures of finite semivariation for families of bounded E^-measurable 
scalar functions /&, by taking the limits of the integrals of n-tuples of 
simple functions uniformly converging to the f^s |Do87] . This inte- 
gral coincides for compact sets with the integral discussed by Morse 
and Transue [Mo56] whose main properties were reviewed in |Ch83 



Among them we must point it out an extension of Lebesgue dominated 
convergence theorem (see example 2.5 and the comments below, Thm. 
2.8 and Cor. 2.9 (iii) in [Ch83]). 

Let Bo(.K'z) denote the Borel a-algebra on the compact space K\. A 
polymeasure 7 on the product of the a-algebras Bo(i^i) x ■ • • x Bo(i^ n ) 
is said to be regular if for any Borel subsets A[ C K h I 7^ k, the set 
function: 

j k (A) = j(Ax, ..., Afe_i, A, A k+1 , . . . , A n ) 
is a signed Radon measure on Ki, I = 1, ... ,n. 
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The space of regular countably additive polymeasures on Bo(Ki) x 
• ■ • x Bo(K n ) will be denoted by rcapm(Bo(i^i), . . . , Bo(K n )) which is 
a Banach space equipped with the semivariation norm. On the other 
hand, we will denote by C n (C(Ki), . . . , C(K n ); R) the space of conti- 
nuous scalar n-linear maps on the Banach spaces C(Ki). Then there ex- 
ists a natural isometric isomorphism between C n (C(Ki), . . . , C(K n ); M), 
and the space of regular countably additive polymeasures rcapm(Bo(i ; fi), 
Bo(AT n )) [Bo98] . 

We will use throughout the rest of this paper a consistent multi-index 
notation. We introduce the symbols V r /z s = V^V^ 2 ■ ■ • V^ n /i Sl ... Sn , 
where V; denotes the difference operator on the Ith component, V;/ik = 
/^k — /^k+i, • Then we obtain easily: 

V>.= >(-D" I , )/':,-! 



1=0 v 7 



where r, s, 1 denote multindexes of length n. We define the Bernstein 
coefficients A(k- m ) of a n-linear functional L as: 

(2.2) 

^ ) t mi (l - tp~ m \ . • • , ( ^ ) t m "(l - t) k "' m ^ = A (k;m) . 

/ i, \ 

7k— m_ 



If L is the functional defined by Uk, then A(k ;m ) = I ) V m /U r 

A function \x on Nq is called completely monotone if V r u > for 
all r. Because of eq. (I2.2p this is equivalent to the positivity of the 
functional L^, hence to the existence of a Radon measure solving the 
classical Hausdorff problem of moments and to the r-positivity of the 
function /j, (see [Be84j for a thorough discussion of these results, Thm. 
4.6.4). We will introduce now two notions of uniform boundedness 
for a multi-index sequence «k that will characterize the solutions of 
the classical and weak multilinear Hausdorff problems in the situation 
where the function \x is not completely monotone. 

Definition 2.1. A multi-index sequence Uk is said to be bounded with 
constant C > if : 



(2.3) ^IV;m)l<C Vk>0, 

m=0 

and Uk is said to be weakly bounded with constant C > if : 



(2.4) 



7 , a m A (k;m) 



m=0 



<C, Vk > 0, 



where aj^ = a^ • • • a^ n , for all a^ ; such that |a^ ; | < 1, Z = 1, . . . , n. 
It is clear from the definitions that condition (12. 3 j) implies (12. 4p . 
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Theorem 2.2. A multi-index sequence fi^, k e Ng solves the weak 
H aus dor ff problem if and only if it is weakly bounded. 

Proof: If /ik is a solution of the weak Hausdorff multilinear moment 
problem, eq. (jl.5|) . then there exists a continuous multilinear functional 
L on C[0, 1] such that L(t k ) = /i k . Let aj^ = a^ 1 • ■ ■ a^ 1 , k > 0, with 



ah,\ < 1, Z 



k 



n. Then: 



7 , a m^ 



(k;m) 



m=0 



< 



Z2 a m L (\k umi )(t), • • • , \k n ,m n )(t)) 
m=0 

/ "T r^n 

L ( E < 1A (fei,mi)(*)>"--> E <C A (*»,™»)(*) 

r. 



vmi=0 



m n =0 



Z=l m;=0 



Moreover, because 

t) mi = 1, we reach the conclussion 



Efc L 
m i= a ki A (ki,mi) 



OLu l \(k, .m,\(f) < Xlm 



m ; =0 



fe 

??!, 



t^- m '(l 



Conversely, if we assume that the multi-index sequence /ik of length 
is weakly bounded, we will show by induction on n that \L(pi, . . . , p n ) \ < 
2 n C||pi|| 00 • • • HpjjIIoo, for any family of polynomials pi, . . . ,p n . 

For n = 1, because of (II. ip the weakly bounded condition for se- 
quences is equivalent to the condition of boundedness. 

We will assume that if V is a (n — ^-multilinear functional associ- 
ated to the bounded (n — 1) -multi-index sequence /x k , with bounding 
constant C, then 



\L'(qi 



qn-i)\<2 n - L C'\\ qi 



\Qn-l\ 



for any family of polynomials qi, ■ ■ ■ ,q n -i- Let L be the multilin- 
ear functional associated to the weakly bounded n-multi-index se- 
quence /ik- We denote by L p the (n — l)-multilinear functional ob- 
tained by fixing the nth argument of L to be the polynomial p, i.e., 
L p (q u . . . , g n _i) = L(gi,...,g n _i,p). 

Notice that if p(t) is a polynomial of degree r then [Sh70] : 



r-l 



B N ( P )(t) = P (t) + S N (t) = P (t) + J2 ^: 



where B^{p) denotes the Nth. Bernstein polynomial of the function p(t) 
and p r i are polynomials of degree less than or equal to r, not depending 
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on N. Denoting by a!L = oh ■ • -a*:™- 1 with \ah I < 1, it is clear that: 



m'=0 



< 



< 



< 



/ , a m'^B N (p)(Ki,mi(tl), • • • , A fen _ 1 , mn _ 1 (t n _i)) 
=0 
k' N 

/ , / , a m 'P(m n /N)L(X klyrni (ti), . . . , A fen _ 1 , mn _ 1 (t n _i), Ajv,m„(*n)) 

n=0 

k 

/ , a m-^(^fci,mi(tl), • • • , Afe n _ 1 , mn _ 1 (t n _i), Afc n , m „(t„)) 



m'=0m n =0 
k 



< 



m=0 



withk= (fei, . . . , fc„_i,iV), m = (toi, . . . ,m ni ,m n ) and a^ = a^ ll (p(m n /N)/\\p\\ 00 ). 
Hence, 



k' 



7 , a m'-^B ]v (p)(Afe 1 , mi (tl), • • • , Afc n _ 1 , m „_ 1 (i„_l)) 



m'=0 



< 



— HpI 



7 , a m^(k;m) 



m=0 



< ClbHc 



because fj,^ is weakly bounded with constant C. Then the (n — 1)- 
multimoment sequence defined by L Bn ^ is weakly bounded with bound 
CUpIIoo, and by the induction hypothesis we obtain: 

\L Bn ( p )(pu . . .,Pn-i)\ < 2 n_1 C||pi|| 00 • • -\\pn-i 



oo IF oo- 



Similarly we consider now S^(t) = Ym=i n^-- If the polynomials 
p r j have the form, p r ,i(t) = Y7j=o a ijV i then by choosing: 



(2.5) 



a = max{|a/j|}, N > 



a(r - l)(r + 1) 



and using the notations above, we will get: 



(2.6) 



m'=0 



< 



7 , a m' L S N (p){^k 1 ,m 1 {tl), ■ ■ ■ , Xk n - U rn n -i(tn-l)) 
=0 
k' r-1 

/ j fl m' / ^(^i,mi(^l)) • • • , ^k n -i,m n -i(t 



< 



n-\)i 



m'=0 
k' 






r r— 1 



< 



2_^ V 2^ 2_^ J-jjL(\ kurrll (ti), . . . , Afc n _ 1 ,m n _ 1 (^n-l),^ 
m'=0 J=0 J=l 

k' / . , r-1 

2J A/7 L ( A fel,mi(^l 



< 



r + 



DO V~^ 



7 , "m' 



m'=0 



i=l 



• • • 5 A fc n _i,m n _i (^n-ljj AjjV'n)) ■ 
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If we denote by e^- the quantity ttw-YmZi W' conditions eq. (12. 5p . 

imply that \a 3 A < 1. We will consider now the numbers at, = for all 
< m n < j. With these definitions the last term in the sequence of 
inequalities fl2.6|) . can be written as: 



r + 



30 V~^ 

I ^ 

3=0 



k' j 



''Xh 



r + 



tE 

3=0 



m'=0 m„=0 

• • • j Afc„_i,m„_i I'-n-lJ) ■^j,m n {*n)) \ 
k 

/ , a m^(^k 1 ,m 1 (h), • • • , Afc n _ 1 , m , l _ 1 (t n -l), ^j,m n {tn)) 



m=0 



with k = (ki, . . . , k n -i, j) and m = (mi, . . . , m n ). Hence finally we 
obtain: 



7 , a m'-^5 Ar (p)(Afc 1 ,m 1 (il),- • •, Afc n _ 1>mn „ 1 (t n _i)) 



m'=0 



< 



< 



r + 



DO V~^ 

i=o 



7 , a mV;i 



m=0 



< 



r + 



^X)C = C||p|| c 



j=0 



and the sequence of multimoments /4„ = L 5jv(p )(A fclimi (ti), . . . , A fcn „ limn _ i ; 
is weakly bounded with constant CHpH^. 

We conclude the argument by using the induction hypothesis and 
computing: 

\L(pi, . . . ,p n )\ = \L Pn (p 1 ,...,p n -i)\ 

< \ L B N {p n )(pi,---,Pn-l)\ + |Lg^(p n )(pi,...,P„_l)| 

< 2 n C||p 1 || 0O ...||p«IU □ 



It is clear that if the multi-index sequence /Xk is nonnegative, i.e., 
/ik are positive or zero real numbers for all k, then the sequence /ik is 
weakly bounded iff is bounded because, 



m=0 



A 



(k;m)| 



E 

m=0 



A(k;m) < 



sup 

i,|om,|<l 



E 

m=0 



A 



(k;m) 



<c. 



Moreover under these circumstances, it is simple to see that the total 
variation of the polymeasure 7 determined by /ik is finite, hence the 
polymeasure 7 determines a Radon measure on [0, 1]" [BoOlJ. Thus we 
have obtained a particular instance of the fact that Radon bimeasures 
on Hausdorff spaces are extensions of Radon measures on the product 
spaces (see |Be84j . Thm. 1.1.10). 
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3. The strong Hausdorff multilinear moment problem: 
orthogonally additive polynomials 

Given a multi-index sequence /ik, k = (k\, ■ ■ ■ , k n ), k[ = 0, 1, . . ., 
Z = 1, . . . , n, we will call it Hankel if /ik+i ; = £ik+i i+1 , where the multi- 
index 1; is defined as (l;)j = 5ij, for all / = 1, . . . , n. 

Theorem 3.1. Let ^ be a multi-index sequence solving the classical 
Hausdorff problem of moments. Then /ik solves the strong Hausdorff 
problem of moments if and only if /ik is Hankel. 

Proof: Consider the n-linear functional defined on the space V of real 
polynomials on / = [0, 1] by the multi-index sequence /ik- Because the 
multi-index sequence /ik is bounded, then L can be extended to C(I) 
(Thm. I2.2p . We shall denote such extension with the same symbol L. 
The homogeneous polynomial Pl determined by L is orthogonally 
additive. To prove it we notice that L(/i, . . . ,g ■ fi, fi+i, ■ ■ ■ , f n ) = 
L(fi, ■■-,fi,9- fi+i, ■■■Jn) for all /i, ...,/»,£ 6 C(I). In fact we can 
construct a sequence of polynomials p mi , q m converging uniformly to f 
and g respectively (I — 1, ... ,n) on I, hence because /ik is Hankel, we 
have: 

*-'\Pm\ i ■ ■ ■ i Qrn'Pmi j Pm; +1 j ■ ■ ■ j Pm„) *-'\Pm\ ■,■■■■, Pmi j Qm'Pmi +1 j ■ ■ ■ i Vm n )i 

and the conclusion follows because of the continuity of L. 

Now suppose we have two disjoint positive functions f,g on C(I), 
|/| A \g\ = 0. We compute: 

P L (f + g) = L(f + g,...J + g) = J2( n r )L(f, n - r J,gZ,9) = 

= L(f, ...,/) + J2 ( n r ) L ^ /. "" r_1 > /, 9, r ~\ gJ-g) 

r=l ^ ' 

+ L(g,...,g) = P L (f)+P L (g), 

because / ■ g = 0. 

Using the representation theorem for orthogonally additive polyno- 
mials on Banach lattices |Be06j and because the n-concavification of the 
Banach lattice E = C(I) coincides with itself, this is C(/)( n ) = C(I) 
(see also |Ca06] and [Pe05] ) . the polynomial Pl defines a bounded lin- 
ear functional T: C(I) — >• M., 

T(P) = P L (f) = L(f } ---J) 

and then, by Riesz theorem, there will exists a Radon measure \x such 
that T(/ n ) = / 7 /(t)"/i(rft). Hence L(t k \ . . . , t fc ") = J I t kl+ ^ +kn ^(dt), 
and L(/ 1; ...,/„) = fj h(t) ■■ ■ f n (t) n(dt). D 
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4. The weak bilinear Hausdorff problem of moments and 
weakly harmonizable stochastic processes 

We will use now the characterization of polymeasures with compact 
support obtained before to provide a description of a class of weakly 
harmonizable processes. Let us consider an stochastic process X t mod- 
elled on a probability space (Q, E, P) where £ is a cr-algebra on the set 
Q, P is a probability measure on Q, and the map X: R — > L 2 (Q,P), 
X t := X(t), is strongly continuous. We denote by C(t,t') = E(X t X' t ) 
the covariance function of X t . If the second order process X t is weakly 
stationary, i.e., there exists an univariate continuous function $ such 
that C(t,t') = $(£' — t), Cramer-Kolmogorov's theorem shows that 
there exists a stochastic measure £ on R with values on L 2 (Q, P) such 
that the process X t is the Fourier transform of £: 



(4.1) X t = / e lts £(ds), 

Jr 

Moreover if A, B are two Borel sets on R, then 
(4-2) (£(A),£(B)) LW) =M^nB). 

Such processes are called (strongly) harmonizable. An important class 
of second order stochastic processes X t that admit generalized har- 
monic representations are the so called weakly harmonizable and they 
satisfy: 

C{t, if) = [ e~ its ® e u ' s ' 7 (ds, ds'), 

JRxR 

where 7 is a positive definite bimeasure, this is: 

r 

(4.3) 7 (A, B) = j(B,A), J2 WM^, Aj) > , 

for all families of complex numbers Oj and Borel sets A, B, A\ on R. 
Then if X t is a weakly harmonizable process then there exists a har- 
monic representation of the form eq. (14. ip for them, where now the or- 
thogonality condition ( I4.2p is replaced by (^(A),^(B)) L 2^i >P ^ = "f(A, B) 
for any Borel sets A, B (see for instance the review |Ra82j . Thm. 3.2.) 
We will consider a complex regular bimeasure 7 of finite semivaria- 
tion with support in [0, 1] x [0, 1]. Let us call such bimeasures Hausdorff. 
We will consider the Fourier-Stieltjes transform of the bimeasure 7: 

its o At' V 



7 (M')= / e~ lts ®e lts j(ds,ds'). 

JRxR 

The function 7 is bounded by 7(0,0) = 7QO, 1], [0, 1]) = /xoo < HtII < 
00 and the extension of Lebesgue's dominated convergence theorem 
mentioned in section 2.1 shows that the function 7 is analytic in the 
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real plane (t, t') with power series expansion given by: 



An-\-m 

(4-4) 7(M')= V (-lr-p-r^,/"/ 



1 

■n±lm 



n\m\ 

n,m>0 



where the coefficients /i nm , n, m > are the moments of the bimeasure 
7. Moreover the sequence of moments /i nm is weakly bounded because 
of Thm. 12.21 Now a simple argument shows that these conditions char- 
acterize completely the Fourier-Stieltjes transform of Hausdorff poly- 
measures. 

A weakly harmonizable second order stochastic process X t such that 
the support of the stochastic measure £ defining it is contained in the 
interval [0, 1] will be called Hausdorff. Notice that in such a case be- 
cause of eq. (I4.ip . the support of the corresponding bimeasure 7 will 
be contained in [0, 1] x [0, 1]. Now if we are given a arbitrary collection 
of complex numbers a±, . . . , a r , and we compute Ymic=i ®i a k'y(ti, tk) for 
a positive definite bimeasure 7, we obtain: 

r „ r r 

J2 Wkifa, t k )= J2 ^ e ~ ltis ® Yl a ^ tkS ' ^ ds > ds "> - °- 

l,k=l JRxR l=1 k=1 

and we conclude that the function 7 is a positive definite kernel. Hence 
the class of analytic positive definite kernels described above are just 
the covariance functions of Hausdorff weakly harmonizable stochastic 
processes. 

Corollary 4.1. A function of two real variables $(£, s) is the covari- 
ance function of a second order weakly harmonizable Hausdorff process 
X t if and only if is an analytic positive definite kernel on R 2 such that 
the multi-index sequence /i nm = d n+m Q/dt n ds m (0, 0), n,m G No is 
weakly bounded. Moreover, the stochastic process X t will be weakly sta- 
tionary if and only if the multimoment sequence /i nm is Hdnkel, i.e., 

Mn+l,m = AVm+1 f 0r al1 n > m - 

Remark 4.2. Notice that the positivity condition can be dispensed 
with as it follows from the previous discussion that an analytic function 
$(s,t) wiill have the form (14. 4 p for a Hausdorff bimeasure 7 iff its 
sequence of moments is weakly bounded. However, unless the positivity 
condition (14.31) is satisfied it is not possible to reconstruct a Hilbert 
space where the stochastic process would be represented. 
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